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We investigate the applicability of Migdal-Eliashberg (ME) theory by revisiting Migdal's 
analysis within the dynamical mean-field theory framework. First, we compute spectral 
functions, the quasi-particle weight, the self energy, renormalised phonon frequency and re- 
sistivity curves of the half-filled Holstein model. We demonstrate how ME theory has a 
phase-transition-like instability at intermediate coupling, and how the Engelsberg-Schrieffer 
(ES) picture is complicated by low-energy excitations from higher order diagrams (demon- 
strating that ES theory is a very weak coupling approach). Through consideration of the 
lowest-order vertex correction, we analyse the applicability of ME theory close to this tran- 
sition. We find a breakdown of the theory in the intermediate coupling adiabatic limit due 
to a divergence in the vertex function. The region of applicability is mapped out, and it 
is found that ME theory is only reliable in the weak coupling adiabatic limit, raising ques- 
tions about the accuracy of recent analyses of cuprate superconductors which do not include 
vertex corrections. 
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I. INTRODUCTION 

The use of Migdal-Eliashberg (ME) theory for the study and analysis of electron-phonon systems 
is widespread. For conventional superconductors such as lead, where the electron-phonon coupling 
is higher than can be treated with BCS theory, ME theory has been extremely successful for 
understanding the superconducting properties. Recently, a lot of researchers have been interested 
in an apparent kink in the electronic dispersion of cuprate supercondutors as determined from angle 
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resolved photo-emission spectroscopy (ARPES), and their analysis / interpretation typically uses 
the related Engelsberg-Schrieffer (ES) result (i.e. just the lowest order Fock diagram) l|, l3, ISJ. The 
ES result divides excitations into long-lived (coherent) low-energy excitations, and rapidly decaying 
high-energy excitations, with a kink at the phonon energy. ES based analysis of ARPES results 
suggests a large electron-phonon coupling in the cuprates, with estimates of the dimensionless 
coupling constant lying between ~ 0.3 and 1.5 (depending on doping) and very large phonon 



frequencies of ~ 50 meV |3[]. There has also been a development of a maximum entropy technique 
for analysing cuprate superconductors which makes use of Eliashberg theory [4]. The analysis 
in Ref. y determines a dimensionless electron-phonon coupling of A > 1, but this is well above 
the A ~ 1 value where one would normally expect perturbation theory to fail. The region of 
intermediate electron-phonon coupling is also relevant to a number of other materials. For example, 
electron paramagnetic resonance measurements on the manganites support a strong "electron- 



5i]. In light of the current experimental situation, 
y, a full re-analysis of the perturbation theory is 



phonon" coupling leading to Jahn-Teller polarons 
and the importance of the conclusions of Refs. 
of high significance. 

The neglect of vertex corrections suggested by Migdal leads to a theory where an infinite set of 
Feynman diagrams may be summed. According to an analysis carried out by Migdal, this theory 
should be valid in the physical regime of electron-phonon problems, where the phonon energy is 
significantly smaller than the intersite hopping ("Migdal's theorem") specifically that the 

vertex corrections are small when Acjo/c/ ^ 1 (1^0 is the phonon frequency and ej the Fermi 
energy) [s^. On the basis of Migdal's analysis, it is often believed that Migdal-Eliashberg theory 
is applicable above A ~ 1 in the adiabatic limit because of the apparently small size of the vertex 
corrections, even though in general perturbative approaches break down (i.e. the functional form 
of the self-energy becomes incorrect) when the coupling constant becomes large. In this article, we 
revisit Migdal's analysis by examining the large-d limit (local approximation). 

The large-d limit has been an effective workshop for determining the validity of approximate 
schemes. The quantum Monte-Carlo (QMC) solution in the large-d limit ^, ol has been compared 
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12] demonstrating that self-consistent 



with a number of different diagrammatic approaches 
second-order perturbation theory worked better than Migdal-Eliashberg theory at intermediate 
phonon frequencies and strong coupling (ujq ~ 0.2t, U > t). The QMC solution is difficult for very 
low phonon frequency and low temperature. QMC self-energies and Green functions are generated 
along the Matsubara axis. As such, it is not easy to make quantitative conclusions about the point 
of breakdown of ME theory as small differences at individual Matsubara frequencies may result in 
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relatively large differences in the spectral functions. Using an alternative approach, Benedetti et 
al. [3] found the formation of more than one extremum in the path-integral formulation at very 
low phonon frequencies and intermediate coupling, leading to a breakdown of ME theory. Work 
on the finite-dimensional Holstein model also shows deficiencies in ME theory. On the basis of a 
comparison with exact diagonalization results, Alexandrov et al. [14] have shown that ME theory 
may break down at intermediate couplings even in the adiabatic limit. 

In this article, we investigate the Holstein model of electron-phonon interactions [l^, which 
employs a number of approximations. The phonon dispersion is flat corresponding to independently 
moving ions and phonon anharmonicity is neglected, resulting in a Hamiltonian, 

<ij>(7 icr i i 

where (cjo-) create (annihilate) electrons at site i with spin a, Xi is the local ion displacement, 
Pi the ion momentum, M the ion mass, t the electron hopping parameter and the chemical 
potential. A mean-field Coulomb pseudopotential may also be included, but in the normal state, 
it is just absorbed into the chemical potential. In the following, t = 0.5 and all energies (including 
temperature) are measured in units of 2t. 

An expression for the effective interaction between electrons can be obtained by performing 
a linear transformation in the phonon variable to remove the electron-phonon term [3]. The 
resulting interaction is retarded, with Fourier components: 

\Muj^) (w2+w^) (t^2+w^) 
where ojg = 2ttsT are the Matsubara frequencies for Bosons. Taking the limit cjq ^ oo, (7 ^ oo, 
while keeping the ratio g/uoQ finite, leads to an attractive Hubbard model [l?! with instantaneous 
attraction of magnitude U = /Mu^ = W\, where W is the half band-width (note that U is 
related to the bipolaron binding energy for the strong coupling two-electron problem). In the 
opposite limit (wq ~^ 0, M 00, keeping Mujq = k finite) the phonon kinetic-energy vanishes, 
leaving only a static variable Xi representing the phonon subsystem in the Hamiltonian, H = 
—t J2<ii>a '^Icr'^ja + J2iai9^i ~ M)"-icr + f Z]j ^f- Thus the problem looks like that of a single electron 



in a disordered potential. The large-c? limit of this model was exactly solved by Millis et al. 



and extended to deal with long range order by Ciuchi et al [19[ |. Thus, the phonon frequency may 



1^ 



be thought of as a parameter for tuning the level of electronic correlation 

In this paper, we investigate the validity of the Migdal-Eliashberg approach within the dynami- 
cal mean-field theory formalism. As we carry out the self-consistency on the real-axis self-energies, 
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we can investigate the behaviour of ME theory over a full range of temperatures and phonon fre- 
quencies. We begin by calculating self-energies, spectral functions, the quasi-particle weight and 
renormalized phonon frequency. To determine the validity of the theory, we calculate an expression 
for the lowest-order correction to the vertex function, and evaluate its magnitude. We find that the 
theory can also break down in the low frequency regime, contrary to the standard interpretation 
of Migdal's analysis. Finally, we calculate resistivity curves in the regime where the lowest-order 
corrections are small. 



II. MIGDAL ELIASHBERG THEORY IN THE LOCAL APPROXIMATION 



In this article, we use the local approximation or dynamical mean-field theory (DMFT) as a 
way of analysing the Migdal-Eliashberg theory. The self-energy of correlated-electron systems is 
momentum independent in limit of large dimensions 211], and approximately momentum indepen- 



dent in 3D. In large-d, lattice models map onto Anderson impurity models with a self-consistent 



hybridisation |22l |. 

The self-consistent DMFT equations may be obtained by rewriting the action for the model to 
be considered in terms of an effective single-site action 3] 

SeS = - f dr f dr' J2 4A^)GoHr - r')Q^(r') + 5int. (3) 
Jo Jo ^ 

Here QQ{iuJn) plays the role of the host Green function in the equivalent impurity model. If one 
assumes that correlations carried via the bath between electrons entering and exiting a single site 
can be neglected (true in the case of large coordination number), the degrees of freedom associated 



with all but one site can be integrated out. This leads to the following self-consistent equation 



gQ^{iuJn) = iuJn + fJ' + G ^{iuJn) - R[G{iuJn)], (4) 

where G{iLOn) is the site-local (impurity) Green function, and is itself a functional of QQ{iuJn)- 
LOn = 2'KT{n + 1/2) are the Fermionic Matsubara frequencies. R[x\ is the reciprocal function of the 
Hilbert transform, defined as R[D{^)] = ^. The Hilbert transform of the non-interacting density 
of states, P(e), is defined as b{i) = deV{£)/{i - e). 

When electrons move on a tight-binding hypercubic lattice, the bare DOS takes the form of 



a Gaussian 



2l|, T>{e) = exp(— e^/2t^)/t-v/27r and there is no simple expression for the reciprocal 



function. Introducing the modified Dyson equation, 

QQ^{iuJn) = G^^{iujn) + S(iw„), (5) 
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(where T,{iuJn) is the electron self-energy), equation ^ may be rewritten as, 

R[G{iu;n)] = iojn + - ^(zWn), (6) 

and inverted to give an expression for the Green function in terms of the self-energy, 

^ f deV{e) _ 

J luJn + A* - - e 

This also allows the approximation to be interpreted as a course graining of the momentum space 
[2^ . To complete the scheme, a form for the electronic self-energy must be calculated in terms of 
the non-interacting Weiss field. This is normally approximate. Then a self-consistent procedure is 
followed: Compute the Green function from equation ([7]), the Weiss field from equation ([5]) and 
then re-calculate the self-energy until convergence is reached. 

The application of ME theory within DMFT corresponds to computing the self-energy and 
phonon propagator from the diagrammatic equations in figure [T]^ a and b). The self-consistent 
solution of such a self-energy corresponds to the summation of all Feynman diagrams which contain 
no vertex corrections. In the low (non-zero) phonon-frequency limit, Migdal's analysis indicates a 
condition, Uujq <C t^, for the neglect of corrections to the vertex function ^. 

The solution of the Dyson equation seen in figure [T][a) results in the full phonon propagator, 

^2 

^^'^^ = (., + ir?)2-4[l + no(u;)]' 
which is then used for the calculation of the electron self-energy. The phonon polarisation bubble 
(dimensionless self-energy) no(iws) is given from perturbation theory, 

no(iws) = -2UT^G{iuJm)G{iuJs +Wm), (9) 

m 

and may bG analytically continued by introducing tlie spectral representation, GiiujYi) — 
/ dxp(x) / {iujn — x) where p{u)) = Im[G(u; + ^^)]/vr and performing the sum over Matsubara fre- 
quencies to give, 

Im[no(u;)] = -2U / dxp{-x)p{uj - x), (10) 
Jo 

at absolute zero, and 

Mn„(.,r)i = -u^M^,jy-±^^^^^. (11) 

at finite T. The full spectral function is used so that all diagrams with no vertex corrections are 
included, consistent with the proper interpretation of Migdal's analysis. 
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The lowest-order skeleton diagram shown in figure [T]^b) , 

E(ia;„,) = — f/T^ G{iu>n — i(^s)D{iu>s), (12) 

s 

may be analytically continued in the same way as the polarisation bubble to give, 

lm[T.{io)] = U r dxp{x)a{LO - x), (13) 
Jo 

with cr{uj) = lm[D(uj + irj)]/'Tr. A more complicated expression applies at finite temperature, 

^ r^, f°° dx p(uj — x)a(x) 

ln.pi.)i = UcoM^)l^^-^^^l±^^. (14) 

At each stage, the Kramers-Kronig relation is used to compute the real parts of the electron and 
phonon self-energies, for instance Re[no(w)] = P / c?xIm[no(x)]/7r(u; — x), where P denotes the 
principal integral. 

After the self-consistent procedure has converged, physical properties are calculated, including 
the quasi-particle weight (inverse effective mass), 



Z = ^ = 1- — 

m* du 



(15) 

LU = 



the effective phonon frequency Q from, 

fi2-a;2(l + Re[no(0)]) =0, (16) 

and the optical conductivity [25l ]. 



Re[a(c^)] = - r deV{e) H dv p{e,v)p{e,v + oj)[f{v) - f{v + uj)], (17) 

^ J —oo J —oo 

where /(x) is the Fermi-Dirac distribution and p{e, v) = Im[l/(z^ + irj — e — S(i/))]/7r (taking the 
limit, Lo duj, the DC conductivity is recovered.) 

Although DMFT is approximate, the formalism can be expected to work in the non-interacting 
limit, where course graining will give the exact non-interacting DOS for the tight-binding model 
(regardless of dimension). In the opposite limit {t 0), the neglect of loops through the host is 
justified, and the formalism should be exact. Such propagation of correlation is also small when 
the coordination number is high (e.g. FCC lattices). 



III. SPECTRAL PROPERTIES 



We have solved the DMFT equations using the self-energies in equations ()13p and (I14p to find 
electron spectral functions. We show the evolution of these functions with coupling for ojq = 0.125 
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in figure [2]|^a) and ujq = 0.5 in figure [2]||b). The spectral functions show two features. There is a 
low energy peak with a width defined by the phonon frequency, and a high energy shoulder. This 
behaviour can be related to the two limits of the Holstein model. At low energy scales (w < loq), 
electrons interact via virtual processes and the behaviour is essentially Hubbard- like (correlated). 
At frequencies greater than ujq, the spectral weight is reduced, phonons may be created, and static 
behaviour emerges (in the static limit, phonons may always be created as no energy is required 
to deform the lattice). The central peak narrows with increased coupling, until a critical value is 



reached. Here, the theory has a Brinkman-Rice-like transition with a diverging effective mass [26l |. 

In figureO we plot the imaginary part of the self-energy. For low coupling strengths (small U and 
energy scales (|a;| < ujq), lm[S{Lj)] is small because electrons cannot create phonons, consistent with 
the Englesberg-Schrieffer analysis). However, as the coupling increases, electrons can be scattered 
by bipolaron resonances and Im[S(a;)] rises sharply on either side of the Fermi-energy. The gap 
in the weight of the self-energy is frequently assumed in the analysis of ARPES data, but as we 
can see here, when the electron-phonon coupling approaches the band width (i.e. A ~ 1), the 
self-energy has a more complicated form, without the simple picture of coherent and incoherent 
quasi-particles. Indeed, it has been demonstrated that in low dimensions, vertex corrections are 
required to achieve a sharp discrimination between coherent and incoherent dressed electrons j^]- 
The onset of this "transition" can be studied by examining the inverse quasi-particle mass 
(quasi-particle weight), shown in figure H] for two values of phonon frequency. As the coupling 
increases, Z becomes smaller and eventually vanishes at a critical value of coupling (the effective 
mass diverges). 

We also examine the phonon spectral function in figure [5] for loq = 0.125 and various couplings. 
As the coupling is increased, the phonon modes soften (figure [6]). The effective frequency does not 
tend to zero as quickly as the quasi-particle weight. We will revisit this point later in this article. 

It is clearly the case that ME theory does not correctly describe the strong coupling limit. It is 



18l | and from approximate "iterated-perturbation 



known from the exact solution of the static limit 
theory" [2g] and QMC [29|] solutions of the Hubbard model, that sub-bands should form at strong 



coupling. As our calculations demonstrate, this sub-band formation is not properly reproduced 
within ME theory and this makes it likely that effects of vertex corrections become significantly 
more important at strong coupling. Since this is at variance with the traditional interpretation of 
Migdal's analysis, we are motivated to re-examine the vertex function in the next section. 
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IV. BREAKDOWN OF MIGDAL ELIASHBERG THEORY 



The neglect of the lowest-order vertex correction shown in figure [D^c) and all higher-order 
corrections is central to Migdal-Eliashberg theory. In this section, we compute the lowest-order 
correction and use this to define the region of validity for ME theory. In this sense, we are 
revisiting Migdal's analysis to understand why there is a contradiction between low frequency 
results computed with ME theory and advanced numerical methods. 

The ratio of first to zeroth order vertices at finite temperature may be written as 



As temperature tends to zero, this sum becomes an integral which may easily be evaluated (x 
and uj are defined as in figure [T|) , 



In order for vertex corrections to be unimportant, the ratio Ti/Tq must be small. As an example, 
we show the lowest order correction to the vertex function for (a) loq = 0.0556 and U = 0.214 (b) 
loq = 0.75 and U = 0.214 in figure [71 For the larger phonon frequencies, the ratio is clearly greater 
than 10% and the approximation is expected to be significantly changed by the inclusion of vertex 
corrections. For the weak coupling U shown here, the corrections are much less pronounced for 
small phonon frequencies. 

In Figure [8] we show the magnitude of the central peak of the vertex correction (always the 
largest part of the function) for a range of couplings at loq = 0.0556. Between U = 0.54 and 
U = 0.55 the vertex correction diverges, and ME theory clearly breaks down. As it is necessary to 
pass through the divergence to reach U > 0.545, the theory is not applicable for strong couplings. 

Figure El shows the lowest value of U' = U/{t + U) for given u'q = LOQ/{t + loq) at which the ratio 
Fi/Fo exceeds 0.1 (the primed quantities are chosen so that U' = 1 represents an unprimed value 
U = oo, and a primed value of U' = 0.5 represents an unprimed value U = t). Below this value ME 
theory can be expected to give accurate results, while the theory begins to break down above this 
line. For ujq ^ 0.2 (loq ^ 0.125) the standard Migdal criterion works well and correctly predicts 
the region in which ME theory is applicable. However for frequencies uj'q ^ 0.2 (uq ^ 0.125), the 
Migdal criterion misses the divergence in the vertex correction associated with the divergence in 
the effective mass. This breakdown happens at f7 ~ t or A ~ 1. The line of breakdown due to this 
vertex divergence levels off as the phonon frequency approaches zero and does not scale as 1/U 




s 



(18) 




(19) 
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(Migdal's criterion). These results imply that there is a breakdown in ME theory in the adiabatic 
limit. 

We note that at high phonon frequencies {U,t <^ ujq), our line does not tend to zero. This is 
because the contribution ri/ro(a; ^ 0,a;' ^ 0) does not tend to infinity. Instead, the function 
gets wider as the phonon frequency increases (see figure [TJ panel b), and some combination of 
the magnitude and width of the vertex is probably a better measure of the importance of vertex 
corrections. At large loq, the breakdown is quite different to that at low phonon frequencies as 
the second order diagrams (one of which contains a vertex correction) develop the same functional 
form and magnitude. However the breakdown at U = 0.1 {W = 0.17) is not too bad an estimate. 
At this value of U, there is only a small mass renormalisation duejto vertex corrections. 

Our results are consistent with that of Benedetti and Zeyher [13], which was obtained with a 
different technique using a semi-circular density of states (Vsci^) = \/4t^"^-^/27rt^), and predicted 
breakdown for Xsc ^ 1-25 in the extreme adiabatic limit. We have included their result in the 
figure, choosing the bandwidth parameter by matching the bare DOS at the Fermi energy for the 
two cases (the conversion factor is A = P5C'(0)Asc'/P(0) = XscVW^) which is reasonable when 
considering low energy excitations at half-filling, and corresponds to a critical A > 0.997 for the 
Gaussian DOS used here (see the diamond in figure[9]). Millis et al. [3] also predicted a breakdown 
in ME theory at C/ ~ t from their solution of the static limit of the infinite-dimensional Holstein 
model, consistent with the breakdown that we have found. The solution of the static limit predicts 
the formation of sub-bands, which is something not reproduced in the ME solution. The lack of 
subbands in ME theory demonstrates that higher order diagrams are essential for the description 
of the sub-band formation, and therefore for the description of the strong coupling limit. 

We suggest that the breakdown of ME theory should be understood in the following heuristic 
manner. As the coupling increases, the host spectral function (^o) narrows and spectral weight 
is shifted towards lower energy scales. There is then an effective bandwidth, teflf < t, for quasi- 
particles close to the Fermi surface. The condition for the applicability of ME theory then becomes 
C/r2 ^ t'^g ^ Z^. As seen in figure [5l there is also a reduction in the effective phonon frequency 
with increasing interaction strength. Although the renormalisation of the phonon frequency helps 
to drive against the transition, we note that the band narrowing effect is much stronger, as can be 
seen by comparing figures S] and [H 
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V. RESISTIVITY CURVES 



Using the finite temperature form of the self-energy, it is possible to calculate resistivity curves. 
These are shown in figure [10] at various electron-phonon coupling and phonon frequency. At high 
temperatures, the tendency is to linear behavior, with a gradient depending on U, but independent 
of uq. The effect of phonon frequency is most dramatically seen at T ~ ujq/2, where a point 
of inflection is seen, before the curve tends to low temperature behavior, consistent with a 
weakly renormalized electron gas. We note that long range order has not been considered in this 
study, so no superconducting transition is seen. The negative curvature of the resistivity curve 
at intermediate temperatures is of interest, since that curvature has recently been interpreted as 
an onset of the Mott limit, but may have a different interpretation as an intermediate coupling 
phenomenon [30]. 



VI. CONCLUSIONS 



We have computed spectral functions, resistivity curves, self-energies, quasi-particle weight and 
effective phonon frequency for the Holstein model at half-filling using Migdal-Eliashberg theory, 
within the DMFT framework, and revisited Migdal's analysis of the strength of the vertex cor- 
rection. By analysing the first order vertex correction, we have defined the region of applicability 
of ME theory. We find that ME theory breaks down at intermediate coupling in the adiabatic 
limit, showing that Migdal-Eliashberg theory should only be trusted at weak coupling within this 
framework. The coupling at which breakdown occurs corresponds to a divergence in the effective 
mass, indicating that both long range order and vertex corrections should be included to correctly 
describe the strong coupling regime. 

The question remains: What was lacking from the analysis of Migdal, which estimated the 
magnitude of the vertex corrections. This can be understood from the magnitude of the kinetic 
energy. In the calculations, the kinetic energy can be seen to be decreasing due to the localisation 
at the non-analyticities. Thus in Migdal's estimate, the renormalised kinetic energy, and not the 
bare kinetic energy should have featured. The renormalised phonon frequency should also have 
been used, but it tends to zero more slowly than the kinetic energy. At the non-analyticity, where 
the kinetic energy tends to zero, it can be seen that a vertex divergence is expected from the 
modified version of Migdal's estimate. 

Our results are timely because of a recent upsurge of interest in the role of electron-phonon 
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coupling in the cuprate superconductors. There are many experimental results which have been 
analysed and interpreted using ES/ME theory, in spite of the large coupling constant that has 
been predicted from that analysis. Many groups are using Migdal's analysis to justify their claims, 
without a careful consideration of the self-consistent effect of coupling on the kinetic energy of the 
electrons (polarons). Exact numerical analysis shows that large electron-phonon couplings signifi- 



cantly reduce the kinetic energy of electrons [3l|] and thus affect the internal consistency of Migdal's 
analysis. Therefore, the any very large A determined from ES style analysis of experimental results 
is not consistent with a theory that neglects vertex corrections (i.e. a more detailed analysis with 
higher order effects included is necessary). Without that analysis, the only conclusion that can be 
reached is that the coupling is large, but no reliable value for that coupling can be determined. The 
main conclusion of a breakdown is not expected to change in ID, 2D or 3D. In three dimensional 
systems, the role of spatial fluctuations should only be significant very close to the bi-polaron 
instability. The role of spatial fluctuations on the spectral properties in 2D has been analysed at 
weak coupling, and is found to lead to quantitative differences in the coherent excitations 27l |. 



Spatial fluctuations have a major role in the superconducting state, and it has also been shown 
that the Eliashberg approach to superconductivity is inadequate for optic al p honon mediated d- 
wave superconductivity 



32( 1 ■ and incomplete for s-wave superconductivity 33|. We therefore urge 
researchers to carefully analyse the internal consistency of their theories when dealing with strong 
coupling materials. 
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FIG. 1: (a) Dyson's equation for the phonon Green function to second order when vertex corrections are 
neglected, (b) Renormalized Fock diagram, (c) The first-order correction to the vertex function, ix is the 
frequency of the emitted phonon and ioj is the frequency of the incoming electron. Neglect of this diagram 
is central to Migdal's theorem. Thick and thin lines represent the full and bare Green functions respectively. 
Wavy lines represent phonons and straight lines electrons. 
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FIG. 2: Spectral functions of the Holstein model computed using Migdal-Eliashberg theory for ujq = 0.125 
(panel a) and luq = 0.5 (panel b) at T = 0. A bi-polaronic resonance forms at zero frequency. No upper and 
lower sub-bands are formed but spectral weight is shifted away from the Fermi-energy. The central peak 
narrows with increasing U, corresponding to a divergence in the effective mass. The general form of the 
curves is similar at all frequencies. 
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FIG. 3: Sclf-cncrgy of the Holstcin model solved using ME theory for luq = 0.125 (panel a) and t^o = 0.5 
(panel b) at T = 0. As the coupling, U — WX (see discussion after eqn 2) increases, Hubbard-likc behaviour 
is seen at low energy scales, and FK-like behaviour at \uj\ > luq ^ 0.2. Note how the gap in the self-energy at 
weak coupling (which is central to the Engelsbcrg SchriefFer approach) fills up as coupling is increased. This 
is due to the increasing importance of higher order diagrams at stronger couplings. This shows that one 
should be careful about interpreting results from the cuprates (A ^ 1) using an ES approach (this type of 
analysis is typically used on ARPES results), since the ES form for the self-energy (i.e. gapped at \oj\ < lvq) 
can only be seen for very small A. 
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FIG. 4: Quasi-particlc weight, Z, computed using ME theory as a function of U (resuhs for phonon fre- 
quencies loq = 0.125 and ojq = 0.5 are shown). Z is strongly renormahsed for intermediate U . 
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FIG. 5: Dependence of the phonon spectral fimction on coupling at ujq = 0.125. The effective phonon 
frequency (location of maximum) and particle lifetime (inverse width of peak) are reduced with increasing 
coupling. When U = 0.5 the curve is skewed in such a way that it is not Lorenzian, and the excitation no 
longer has the properties of a single phonon. 
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FIG. 6: Coupling dependence of the renormalised phonon frequency defined by — llJq{1 + Re[no(il)]) — 0. 
The renormalisation of the effective phonon frequency is not as strong as that of the electronic inverse 
mass (quasi-particle weight Z) shown in figure IH which goes to zero at the transition. The inset shows the 
quasiparticle lifetime given by -y/Im[no(il)]. When the value of A tends to 1, the phonons can no longer be 
treated as single particle excitations. 




FIG. 7: The first order correction to the vertex function when (a) ujq = 0.0556 and U — 0.214 (b) ujq = 0.75 
and U — 0.214. The central maximum in (b) shows a ratio of over 10% and one should expect significant 
corrections to quantities computed within ME theory. 
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FIG. 8: Divergence in the vertex correction. At the divergence Migdal's analysis breaks down. Points for 
larger values of U (grey line) are shown in order to highlight the divergence but the ME treatment used 
here is not valid above the divergence, since the non-analyticity marks the limit of the expansion in A. 




(Oo/(t+G)o) 

FIG. 9: Validity of ME theory as a function ofU' = U/{U + t) and u'q = ljo/{u)o + 1), chosen so that a value 
of 1 represents U = oo, etc. Also shown is the result of Benedetti and Zeyher for very low phonon frequency 
(filled diamond). Above the line, the ratio of the first order vertex correction to the bare vertex, Ti/Tq, 
exceeds 10% and ME theory is no longer strictly valid. At low frequencies (adiabatic limit), breakdown occurs 
at smaller coupling strength than expected as a result of the divergence found in the vertex corrections. This 
is in contrast to the interpretation of Migdal's analysis which is often used when low frequency experimental 
data is analysed, and is the main result of this article. 
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FIG. 10: Resistivity curves for various couplings and temperatures. The resistivity is in units of e^V/ha^ 
with V the unit cell volume and a the lattice cell spacing. The temperature at which the shoulder appears is 
defined by the phonon frequency, and the amplitude is defined by the coupling. Low temperature behaviour 
is according to a weakly renormalized electron gas. High temperature behaviour is linear. Note the 
negative curvature at intermediate temperature. 



